Gemini and Chandra observations of Abell 586, a relaxed 

strong-lensing cluster. 
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ABSTRACT 

We analyze the mass content of the massive strong-lensing cluster Abell 586 (z = 0.17). We 
use optical data (imaging and spectroscopy) obtained with the Gemini Multi-Object Spectrograph 
(GMOS) mounted on the 8-m Gemini-North telescope, together with publicly available X-ray data 
taken with the Chandra space telescope. Employing different techniques - velocity distribution 
of galaxies, weak gravitational lensing, and X-ray spatially resolved spectroscopy - we derive 
mass and velocity dispersion estimates from each of them. All estimates agree well with each 
other, within a 68% confidence level, indicating a velocity dispersion of 1000 - 1250 km s _1 . The 
projected mass distributions obtained through weak-lcnsing and X-ray emission are strikingly 
similar, having nearly circular geometry. We suggest that Abell 586 is probably a truly relaxed 
cluster, whose last major merger occurred more than ~ 4 Gyr ago. 



Subject headings: galaxies: clusters: individual: Abell 586 
lensing - cosmology: observations - dark matter 



X-rays: galaxies: clusters - gravitational 



1. Introduction 

In a Universe where structures are formed hi- 
erarchically, as the current ACDM paradigm pre- 
dicts, larger objects are formed through merging 
and/or accretion of smaller systems and the most 
massive structures should be the youngest. In 
the present era of evolution of the universe, mas- 
sive clusters of galaxies occupy the upper limit 
of the mass function of (nearly) virialized struc- 
tures. This special characteristic makes clusters 
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privileged objects for the study of structure for- 
mation and the nature of dark matter, which is 
believed to play a major role in this process (e.g. 
White & Rees 1978; Kauffmann et al. 1999). 

Cluster masses can be measured by several 
techniques, each one relying on different princi- 
ples and simplifying assumptions, Consequently, 
their biases are also different. Two techniques 
the kinematics of the member galaxies and the X- 
ray emission from the hot gas that fills the intra- 
cluster medium (ICM) - assume dynamical equi- 
librium. On the other hand, gravitational lens- 
ing, in its strong and weak regime, does not re- 
quire such an assumption, depending directly of 
the cluster surface mass (see the reviews by, e.g., 
Fort & Mcllier 1994; Mcllier 1999). 

The use of different techniques to measure mass 
distributions will likely give different mass esti- 
mates. This was first demonstrated by Miralda- 
Escude & Babul (1995), who have found a system- 
atic difference by a factor of ~ 2-3 between strong 
lensing and X-ray mass measurements, with the 
second method systematically producing smaller 
values. Allen (1998) claimed to have solved this 
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problem after comparing results for clusters with 
and without cooling flows. He found that, for the 
former, X-ray and lensing mass estimates tend to 
agree, whereas for the latter the opposite happens. 
With the assumption that cooling flow clusters are 
dynamically more evolved, he concluded that this 
discrepancy is due to non-thermal processes (such 
as merger of sub-components) affecting the ICM 
in the central regions of non-cooling flow clusters. 

Gravitational lensing is often claimed to be a 
more reliable method for cluster mass estimation, 
because it does not depend on equilibrium as- 
sumptions about the dynamical state of the clus- 
ter. Among dynamical mass estimators, the X- 
ray emission of the intra-cluster gas is preferred 
to galaxy dynamics, since the gas is a collisional 
component that reaches equilibrium faster than 
galaxies. In other words, gas is considered a better 
tracer of the gravitational potential than galaxies 
(see, e.g., Sarazin 1988, Sect. 5.5). 

On the other side, there are growing evidences 
that lensing masses for some clusters can be over- 
estimated due to the presence of other massive 
components along the line of sight (Metzeler et al. 
1999; Hoekstra 2003). This seems to be the case 
for the strong lensing clusters A2744 (AC118) (Gi- 
rardi & Mezzeti 2001; Cypriano et al. 2004) and 
CL0024+16 (Czoske et al. 2002; Kneib et al. 2003). 
In these cases, distance information is needed (e.g., 
from redshift surveys) to disentangle clearly two or 
more components along the line of sight. 

Therefore, the use of a multi-technique ap- 
proach to investigate galaxy clusters (e.g., Smail 
et al. 1997; Valtchanov et al. 2002; Ferrari et al. 
2003; Proust et al. 2003; Smith et al. 2003, 2004) 
is ideal for the study of these objects. Not only it 
does allow more reliable mass estimations, but it 
may also offer several hints on the physical state 
of such complex and evolving systems. The goal of 
this paper is to apply this approach to investigate 
the mass distribution and the dynamical state of 
the cluster Abell 586. 

Abcll 586 is a Bautz-Morgan type I cluster, at 
z = 0.17, with Abell richness class 3. It presents 
copious X-ray emission (Lx = 1-11 x lO 45 /^ 1 
erg s _1 ; Ebeling et al. 1998), with its peak co- 
incident with the position of the brightest clus- 
ter galaxy (BCG). We have chosen this particular 
cluster for this study because of its high X-ray 
luminosity and the presence of strong lensing fea- 



tures, both implying a high mass content. 

Dahle et al. (2002) reported a long faint blue 
arc at the northwest of the cluster central galaxy, 
and our optical observations with Gemini revealed 
a fainter arc at the opposite direction with respect 
to the BCG. 

Recent X-ray studies of this cluster (Allen 
2000; White 2000), using ROSAT images and 
ASCA spectra, have found ICM temperatures 
ranging from 6.1 to 8.7 kcV, depending on the 
adopted model. A dc-projection analysis made 
by Allen (2000) predicts a velocity dispersion 
of 1050i25okm s _1 for this cluster. However, a 
weak-lensing mass measurement by Dahle et al. 
(2002) leads to a much larger velocity dispersion: 
a = 1680^ km s -1 . 

In this paper we present new optical and X-ray 
observations (Section 2) that allowed us to esti- 
mate the cluster mass based on galaxy dynamics, 
weak and strong lensing, and ICM temperature 
and surface brightness profiles (Section 3). In Sec- 
tion 4 we compare and discuss the results obtained 
with the different methods and, finally, we summa- 
rize our results and present our main conclusions 
in Section 5. We adopt hereafter H n = 70 /170 
km s- 1 Mpc -1 , Q M = 0.3 and Q A = 0.7. At the 
distance of Abell 586, one arcsecond corresponds 
to 2.9 kpc. All quoted uncertainties are for a confi- 
dence level of 1-cr, or 68%, unless stated otherwise. 

2. Observations and data reduction 

In this section we present the optical and X-ray 
data used in our analysis and describe the main 
steps in their reduction. 

2.1. Gemini-North optical data 

All optical observations discussed here were 
obtained using the Frederick C. Gillett Telescope 
(Gemini North) at Mauna Kea operating in queue 
mode. Imaging and multi-object spectroscopy 
were carried out with GMOS (Hook et al. 2002). 
Image and spectroscopic basic reductions (de- 
biasing, flat-fielding, wavelength calibration, etc.) 
were done in a standard way, with the GMOS 
package running under the IRAF environment. 
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2.1.1. Imaging 

We observed the cluster Abell 586 on two oc- 
casions. The first time (period 2001B) as part of 
a survey for gravitational arcs in 8 clusters with 
X-ray luminosities larger than 10 45 /i^ 1 erg s" 1 in 
the BCS catalog (Ebcling ct al. 1998). 

This imaging consisted of three single exposures 
with the g', r' and i' Sloan filters (Fukugita et al. 
1996), the integration times being 300, 250 and 
250 seconds, respectively. Atmospheric conditions 
were nearly photometric and the seeing, as mea- 
sured by the FWHM of point sources, was ~0.7". 
With these observations we detected gravitational 
arcs in Abell 586, that had already been reported 
by Dahle ct al. (2002). 

During Gemini period 2002B, follow-up obser- 
vations of the cluster were performed, compris- 
ing deeper r' imaging and also multi-object spec- 
troscopy. The total exposure time for imaging was 
20 minutes (4 x 300s). These new images were 
taken in photometric conditions, and the seeing of 
the combined image is again ~0.7". 

In both cases we have used GMOS with 2x2 
binned pixels, leading to a pixel size of 0.145". 
The observed field of view has 5.5' on a side. 

We use the program SExtractor (Bcrtin & 
Arnouts 1996) to build galaxy catalogs, adopting, 
as detection criterion, that objects should have 
at least 10 contiguous pixels with values above 
the background plus 1.5 times its dispersion. For 
the first run, the galaxy catalogs start to become 
incomplete (i.e., the logarithmic number counts 
start to departure from a linear behavior) at 22.5, 
22.5 and 22.0 mag, for g', r' and i' bands, respec- 
tively. The magnitude completeness limit for the 
second run data is r' = 24.5 magnitudes. 

2.1.2. Spectroscopy 

The GMOS spectroscopic observations were 
done using a 400 line mm -1 red optimized grating, 
with a central wavelength of 7000A. This configu- 
ration resulted in spectra with a resolution of 7A, 
or nearly five spectral pixels (as measured from 
the FWHM of calibration lamp lines), covering 
the range 5000-9200A. 

The observations were done with two masks, 
each one with roughly 30 slits, with single integra- 
tion times of 34 minutes each. Slits in the central 



part of the mask were placed over the large grav- 
itational arcs and other candidate lensed objects. 
The remaining slits were positioned over bright 
galaxies in the GMOS field. 

Rcdshifts for galaxies with absorption lines 
were determined using the cross-correlation tech- 
nique (Tonry & Davis 1979) as implemented by 
the software RVSAO (Kurtz & Mink 1998). We 
have used as templates the spectra of the NGC 
galaxies 1700, 1426, 3096, 4087, 4472, and 4751, 
observed with the 2.5-m Du Pont telescope of 
the Las Campanas Observatory, and a synthetic 
spectrum built from the stellar library of Jacoby, 
Hunter & Christian (1984). 

All these templates were observed and reduced 
by Rodrigo Carrasco. More details can be found 
in Carrasco, Mendes de Oliveira & Infante (2005). 
Only object-template matches with the correlation 
coefficient parameter R > 3 were considered suc- 
cessful. Typical formal errors in the radial veloci- 
ties are smaller than 30 km s _1 . 

For galaxies having prominent emission lines, 
the redshifts were also determined directly from 
these lines, but these measurements were adopted 
only when the velocity determinations by cross- 
correlation were doubtful (R < 5). In all cases, the 
differences between the two determinations were 
both results are smaller than 150 km s _1 . In total, 
we obtained 44 reliable redshifts. The resulting 
redshifts for the cluster galaxies are summarized 
in Table 1, and for the other galaxies in Table 2. 

Objects are identified using their J2000.0 coor- 
dinates, in the same fashion as Cohen & Kncib 
(2002), so that C088.3712 has coordinates R.A. 
= 07 h 32 m 08. s 8, DEC. = +31°37'12". The prefix 
"C" in the name means that the galaxy is a cluster 
member, otherwise the prefix is a "G" . 

2.2. Chandra X-ray observation and data 
reduction 

Abell 586 was observed in September 2000 by 
the Chandra satellite in a single 11.83 ksec expo- 
sure with ACIS-I camera (observation ID 530, P.I. 
Leon Van Speybroeck). The data were taken in 
Very Faint mode with a time resolution of 3.24 sec. 
The CCD temperature was — 120°C. The data was 
reduced using CIAO version 3.0.1 1 following the 

http: //asc .harvard . edu/ciao/ 
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Table 1 

Spectral data for Abell 586 galaxies. 



( 1 \ 

I 1 ) 
Name 


(2) 

RA (2000) 


(3) 

DEC (2000) 




(4) 

z 


(5) 
R 


(6) 

r (Ar5 Mag.) 






(') 

Emission lines 




C088.3712 


07 


32 


08.81 


31 


37 


12 


5 


o 


1686 


7 


7 


18.82 


H/3, 




fOIl [NIL 1 Ha 


[Sill 


C093_3835 


07 


32 


09.33 


31 


38 


35 


1 





1687 


6 


7 


20.12 










C094.3542 


07 


32 


09.39 


31 


35 


42 


1 





1739 


6 


2 


19.75 










C107_3616 


07 


32 


10.68 


31 


37 


24 


8 


o 


1744 


13 


7 


18.84 










C107.3725 


07 


32 


10.74 


31 


36 


16 


5 





1728 


7 


4 


18.86 


H/3, 


[OIIIl. 

L J ? 


[Oil, [NIL]. Ha, 

L J" L J 5 7 


[Sill 

J 


C118_3731 


07 


32 


11.79 


31 


37 


31 


2 


o 


1699 


12 


9 


18.34 










C123_3605 


07 


32 


12.35 


31 


36 


04 


7 


o 


1665 


6 


5 


20.62 










C127.3828 


07 


32 


12.74 


31 


38 


27 


8 


o 


1655 


G 





20.86 










C131.3619 a 


07 


32 


13.13 


31 


36 


18 


6 


o 


1597 






20.14 


H/3, 


[OIIIl 

1.1 ; 


[Hell [Oil [Nil 

L ^ J j L J 5 L 


, Ha, [SII] 


C 137.3845 


07 


32 


13.70 


31 


38 


44 


5 


o 


1714 


12 


8 


18.87 










C144.3717 


07 


32 


14.42 


31 


37 


17 


3 


o 


1730 


13 





18.04 










C156_3838 


07 


32 


15.56 


31 


38 


37 


7 


o 


1767 


10 


4 


19.65 










C157.3723 


07 


32 


15.65 


31 


37 


22 


7 


o 


1802 


14 


4 


18.53 










C163.3727 


07 


32 


16.32 


31 


37 


27 


2 


o 


1698 


10 


8 


19.75 










C169_3839 


07 


32 


16.86 


31 


38 


39 


3 


o 


1683 


14 


4 


18.54 










C171.3650 b 


07 


32 


17.07 


31 


36 


50 





o 


1738 


5 


9 


18.23 


H/3, 


[OIIIl 


[OI], [Nil], Ha, 


[Sill 


C176.3744 


07 


32 


17.56 


31 


37 


43 


8 





1658 


5 


5 


20.41 










C177.3857 


07 


32 


17.69 


31 


38 


57 


3 





1679 


13 





18.06 










C186.3846 a 


07 


32 


18.64 


31 


38 


46 


5 





1673 






20.96 c 


H/3, 


[OIII], 


[Nil], Ha, [SII] 




C187.3846 


07 


32 


18.71 


31 


38 


4G 








1595 


4 


9 


20.96 c 










C197.3726 


07 


32 


19.74 


31 


37 


26 


3 





1754 


7 


2 


21.70 










C223.3818 


07 


32 


22.27 


31 


38 


17 


9 





1720 


16 


3 


19.02 










C231.3752 


07 


32 


23.08 


31 


37 


51 


9 





1700 


13 


1 


18.59 










C233.3800 


07 


32 


23.26 


31 


38 


00 


3 





1780 


6 


2 


20.68 










C257.3708 


07 


32 


25.75 


31 


37 


07 


7 





1758 


8 


8 


19.61 










C266.3652 


07 


32 


26.61 


31 


36 


52 


2 





1666 


4 


4 


21.38 










C268.3631 


07 


32 


26.77 


31 


36 


31 








1714 


14 





20.03 










C273.3618 


07 


32 


27.29 


31 


36 


17 


9 





1715 


8 





19.58 










C276.3652 


07 


32 


27.57 


31 


36 


51 


5 





1725 


21 


3 


18.52 










C281.3813 


07 


32 


28.14 


31 


38 


12 


6 





1686 


8 


8 


17.67 










C291.3557 a 


07 


32 


29.11 


31 


35 


56 


5 





1710 






20.52 


[OIII], Ha, 


[SII] 





a Rcdshift measured from emission lines. 

b This galaxy presents Ha in emission with narrow and broad (> 2000 km s _1 ) components, being probably a type 1 
Seyfert galaxy. This object is also the source of the point-like X-ray emission seen near the Southwest corner of Figure 
5. 

c The images of these galaxies overlap, and the magnitude in the table corresponds to the sum of both images. 
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Table 2 

Spectral data for non-cluster members in the field of Abell 586. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


. { . 7) 


Name 


RA (2000) 


DEC (2000) 


z 


R 


r' (AB Mag.) 


Emission lines 



G137.3635 


07 32 13.67 


31 


36 35.2 


0.2203 


14.1 


19.64 




G147.3845 


07 32 14.72 


31 


38 45.1 


0.2123 


7.5 


18.78 


[Nil], Ha 


G181.3850 


07 32 18.07 


31 


38 49.9 


0.3050 


5.3 


18.66 




G194.3816 a - c 


07 32 19.40 


31 


38 15.8 


0.6093 




18.86 


[Oil], H/3, [OIII] 


G199.3734 a -^ 


c 07 32 19.92 


31 


37 33.6 


1.4302 




21.87 


[Oil] 


G208.3842 a < b 


07 32 20.76 


31 


38 41.7 


0.8973 




21.89 


[Oil] 


G243.3809 


07 32 24.29 


31 


38 08.7 


0.2453 


3.6 


22.25 




G257.3837 


07 32 25.69 


31 


38 37.1 


0.3050 


7.4 


18.54 




G283.3604 a < b 


07 32 28.33 


31 


36 04.2 


0.8563 




21.06 


[Oil] 


G288.3703 a 


07 32 28.78 


31 


37 03.3 


0.1921 




20.56 


[OIII], [Nil], Ha 


G298.3552 


07 32 29.79 


31 


36 23.3 


0.1912 


4.3 


19.21 


[Nil], Ha 


G298.3623 


07 32 29.80 


31 


35 51.6 


0.2141 


8.1 


19.31 




G305.3734 a 


07 32 30.52 


31 


37 34.1 


0.0596 




19.16 


H/3, [OIII], [Nil], 



a Redshift measured from emission lines. 

b Dcspite the presence of only one emission line, [Oil], the redshift of these galaxies has been confirmed 
by the presence of strong absorption lines, namely Call H and K, Mgll and Fell. 

c Strong lensing features. 
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Standard Data Processing, producing new level 1 
and 2 event files. 

The level 2 event file was further filtered, keep- 
ing only ASCA grades 2 0, 2, 3, 4 and 6. We 
checked that no afterglow was present and applied 
the Good Time Intervals (GTI) supplied by the 
pipeline. No background flares were observed and 
the total livetime is 10.04 ksec. 

We have used the CTI-corrected ACIS back- 
ground event files ("blank-sky"), produced by the 
ACIS calibration team 3 , available from the cal- 
ibration data base (CALDB). The background 
events were filtered, keeping the same grades as 
the source events, and then were reprojected to 
match the sky coordinates of Abell 586 ACIS ob- 
servation. 

We restricted our analysis to the range [0.3- 
8.0 keV], because above 8.0 keV, the X-ray emis- 
sion is largely background-dominated. 

2.2.1. X-ray Imaging 

We have constructed an adaptively smoothed 
image in the [0.3-8.0 keV] band using the CSMOOTH 
tool from CI AO. The exposure map was generated 
by the script merge_all, where we have calcu- 
lated the spectral weights, needed for the instru- 
ment map, using the cluster total spectrum, i.e., 
the spectrum obtained inside a circle concentric 
with the cluster, with a radius of 5 arcmin (about 
the CMOS field of view). We first smoothed the 
raw image, then the exposure map, using the scale 
map produced by the smoothing of the raw image, 
and finally divided the smoothed raw image by the 
smoothed exposure map. The result is shown in 
Figure 1. 

The X-ray cluster emission can be detected up 
to 5 arcmin (0.87h^Q Mpc) and is fairly symmet- 
ric. There is a strong point source at 78 arcsec 
from the center toward the Southwest that we 
identify with an active galaxy C171_3650 (cf. Ta- 
ble 1). 

The absence of significant substructure in the 

2 The grade of an event is a code that identifies which pix- 
els, within the three pixel-by-three pixel island centered 
on the local charge maximum, are above certain amplitude 
thresholds. The so-called ASCA grades, in the absence of 
pilcup, appear to optimize the signal-to-background ratio, 
http : //cxc. harvard . edu/ 

3 http: //exc .harvard . edu/cal/Acis/WWWacis_cal .html 
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Fig. 1. — Chandra ACIS-I X-ray adaptively 
smoothed image of Abell 586 in the [0.3-8.0 keV] 
band, corrected by the exposure map and binned 
by a factor of 4. The center of the cluster is on the 
13 chip (bottom right). There are also a number 
of point sources and a faint, extended source on 
the 10 chip (top left). The square (5.5' on a side) 
centered on the cluster shows the GMOS field of 
view. 

X-ray map suggests that the last major merger 
with clusters with masses larger than 1/4 the 
Abell 586 mass occurred on a timescale longer 
than the cluster relaxation time. Roughly, the 
relaxation time is of the order of the dynamical 
time, Td ~ (Gp) -1 / 2 where ~p ~ 340f2j\/Pc inside 
the virial radius for a ACDM cosmology. Thus, 
Td *~ Gyr and we suggest that the last ma- 
jor merger was at least 4 Gyr ago. Such a rough 
estimate usually agrees with ./V-body simulation 
results (e.g. Roettiger et al. 1998; Rowley et al. 
2004). 

2.2.2. X-ray Spectroscopy 

For the spectral analysis, we have computed the 
weighted redistribution and ancillary files (RMF 
and ARF) using the tasks mkrmf and MKWARF 
from CIAO. These tasks take into account the ex- 
tended nature of the X-ray emission. Background 
spectra were constructed from the blank-sky event 
files and were extracted at the same regions (in de- 
tector coordinates) as the source spectra that we 
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want to fit. 

The spectral fits were done using xspec vll.3, 
restricting to the range [0.3-8.0 keV]. The X-ray 
spectrum of each extraction region was modeled 
as being produced by a single temperature plasma 
and we employed the MEKAL (Kaastra & Mewe 
1993; Liedahl et al. 1995) model. The photoelec- 
tric absorption - mainly due to neutral hydrogen 
- was computed using the cross-sections given by 
Balucinska-Church & McCammon (1992), avail- 
able in xspec. 

The overall spectrum was extracted within a 
circular region of 1.4 arcmin (243 kpc) centered on 
the X-ray emission peak. It was re-binned with the 
GRPPHA task, so that there are at least 30 counts 
per energy bin. This radius was chosen because 
almost all emission is in this region, we can avoid 
the CCDs gaps, and we have all the spectrum ex- 
tracted in a single ACIS-I CCD. Figure 2 shows 
the overall spectrum fitted to an absorbed mekal 
plasma. 



> 



free /V H 








' — — |J - Y || jj lt— y u -y - i , = 


fixed iV H = 5. 15 10 20 cm" 2 









0.5 1 2 5 

channel energy [keV] 



Fig. 2. — Chandra ACIS-I X-ray spectrum ex- 
tracted in the central 1.4 arcmin (245/ifo 1 kpc) su- 
perposed with a mekal plasma spectrum. Top: 
fit with Ah free. Bottom: fit with Ah fixed at the 
galactic value. Below each spectrum the residuals 
are shown expressed as the x 2 contribution of each 
binned energy channel. 

A single-temperature gas fits well the overall 
spectrum; the reduced \ 2 is x 2 /d.o.f.= 132.7/133. 
The temperature thus obtained is kT = 7.8^Q g 



keV (90% confidence level), in agreement with pre- 
vious estimates; the metal abundance (hereafter, 
metallicity) is Z = 0.51^o'i7 ^0 an d * ne hydro- 
gen column density is Ah = 9.llj'3 x 10 20 cm~ 2 . 
The metallicity is higher than the average value of 
1/3-Z© for clusters of galaxies. 

The value of Ah obtained from the spectral fit 
is above the galactic value (Dickey & Lockman 
1990) which is A H = 5.15 x 10 20 cm" 2 for the re- 
gion of Abell 586. Not considering the uncertain- 
ties of the galactic Ah measurements, this excess 
is significant at a 3er level. Such a level of signif- 
icance may indicate that we are indeed detecting 
an X-ray absorption excess. The origin of such an 
absorption is controversial (Allen 2000, and ref- 
erences therein), the most promising hypotheses 
being very cool molecular gas and/or dust grains 
that survives in the intracluster medium (ICM). 

If we fix Ah at the galactic value, we obtain 
a much worst fit, x 2 /d-o.f.= 157.6/134 and a 
higher temperature, kT = 9.8+*;QkcV [this anti- 
correlation between the temperature and hydro- 
gen column density is well-known, cf. Pislar et 
al. (1997)]. As can be seen in the residual plot of 
figure 2, this higher value of x 2 is mainly due to 
a poor fit at low energies (E < 0.7keV). At high 
energies (E > 6.5 keV), the fit is also a bit poorer 
than when Ah is considered as a free parameter. 

2.2.3. Flux and luminosity 

Table 3 summarizes the non-absorbed fluxes 
and luminosities obtained inside the central 1.4 ar- 
cmin field. We give the results in "soft" and 
"hard" bands, as well as the rest-frame bolometric 
luminosity (computed by extrapolating the plasma 
cmissivity from 10 eV to 100 keV). 

Using the empirical relation Lx-Tx obtained 
by Xue & Wu (2000), the bolometric X-ray lu- 
minosity (converted to Ho = 50 km s _1 Mpc -1 ) 
implies kT = 6.7±0.7, slightly cooler but in agree- 
ment (within the error bars) with the spectroscop- 
ically determined temperature. 

3. Data analysis 

3.1. Dynamics of galaxies 

We obtained accurate rcdshifts for forty-four 
galaxies in the GMOS field of Abell 586, whose 
properties are summarized in Tables 1 and 2. 
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Table 3: Non-absorbed fluxes and luminosities in different energy bands. The energy band limits are given 
in keV. 



flux 


flux 


flux 


L x 


L x 


L x 


[0.5-2.0] 


[2.0-10.0] 


bolom. 


[0.5-2.0] 


[2.0-10.0] 


bolom. 


2.8 


5.4 


11.3 


2.08 


4.38 


9.18 



Note: flux units are 10 12 erg cm 2 s 
Luminosity units arc 10 44 erg s~ 



Thirty-one of them, with redshifts between 0.16 
and 0.18, constitute our spectroscopic sample of 
cluster galaxies. In Figure 3 we show examples of 
these spectra, and in Figure 4 their velocity dis- 
tribution. The redshifts of the remaining galaxies 
have either z < 0.07 or z > 0.19. 



5000 5500 6000 




5500 6000 6500 7000 7500 



Observed W 



5000 5500 6000 




5500 6000 6500 7000 7500 



Observed W 

Fig. 3. — Examples of spectra of cluster member 
galaxies. The upper and lower panels shows, re- 
spectively, the spectra of the galaxy with the best 
correlation coefficient (C276.3652) and the worst 
one (C266_3652). Both spectra have not been flux 
calibrated and have been smoothed using a boxcar 
filter with 6.8 A (5 pixels) of length, for the sake 
of clarity. 

We estimated the systemic redshift and the ve- 
locity dispersion of the cluster using the bi- weight 
estimators of the ROSTAT program (Beers, Flynn 
& Gebhartd 1990). The resulting mean redshift is 



Av (10 3 km s- 1 ) 
-3-2-10 1 2 3 
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Fig. 4. — Redshift distribution for the 31 galax- 
ies in Table 1 with 0.16 < z < 0.18. The bin 
size is 500 km s _1 . The tick-marks on the top 
of the panel represent individual galaxy veloci- 
ties. The continuous line is a Gaussian function 
with bi-weighted location and scale values equal 
to (z) = 0.1708 and a = 1161 km s" 1 , respec- 
tively. 

(z) = 0.1708 ± 0.0006 and the velocity dispersion 
is cr = 1161 ± 196 km s" 1 . 

This value of (z) compares very well with the 
value obtained by Bottini (2001) using only 7 ra- 
dial velocities of galaxies with very small projected 
distances to the BCG, 0.170 ± 0.001. His deter- 
mined velocity dispersion, 313 ± 70 km s _1 , is 
much smaller than our value; this is not unex- 
pected given his small sample. 

There is a noticeable gap in the velocity his- 
togram of ~ 1500 km s _1 (in the cluster rest 
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frame) between the two galaxies at z ~ 0.160 
and the others. Excluding these two galaxies 
from our sample, the inferred values are (z) = 
0.1711 ± 0.0006 and a = 977 ± 130 km s" 1 . How- 
ever, we consider that these galaxies should be 
kept in our sample because of the following rea- 
sons. Firstly, these galaxies are not excluded by 
a 3-sigma clipping selection. Secondly, the veloc- 
ity histogram becomes too asymmetric after ex- 
clusion of the two galaxies (the ROSTAT asym- 
metry index increases from 0.48 to 0.96). Finally, 
because the difference between the BCG velocity 
and the systemic velocity which is normally small 
in cD clusters (Quintana & Lawrie 1982), becomes 
larger, increasing from 121 to 221 km s _1 in the 
cluster rest-frame. 

We also verified the effect on the velocity dis- 
persion caused by the 6 emission line galaxies, con- 
cluding that their removal from the calculation 
does not cause any appreciable change in a. 

Despite the internal robustness of the velocity 
dispersion found here, a word of caution is war- 
ranted. All galaxies used in this analysis lie within 
a radius of 570 kpc from the cluster center 
(~ 0.4i? vir ; see Section 3.3). In order to under- 
stand more clearly the dynamics of this cluster, a 
larger set of velocities obtained over a larger field 
is required, as has been clearly demonstrated by, 
e.g., Czoske et al. (2002). 

3.2. Weak Gravitational Lensing 

3.2.1. Galaxy shape measurements 

To estimate the cluster mass distribution 
through weak lensing, it is necessary to measure 
accurately the cllipticity of the background galax- 
ies, which includes the effects of distortions intro- 
duced by the cluster gravitational potential. 

The determination of galaxy shapes was per- 
formed using the method described in detail by 
Cypriano et al. (2004). Here we only summarize 
the main steps of this procedure. 

The program im2shape (Bridle et al. 2002) is 
the basic tool we used for shape measurements. 
This program uses Baycsian methods to fit as- 
tronomical images as a sum of two-dimensional 
Gaussian functions with elliptical bases. Each of 
them is fully defined by six parameters: two Carte- 
sian center coordinates, ellipticity, position angle, 
semi-major axis size and the amplitude. More- 



over, IM2SHAPE deconvolves the fitted result by 
using a PSF, given also in terms of a sum of two- 
dimensional Gaussians, thus recovering an unbi- 
ased and accurate estimation of the object's shape. 

IM2SHAPE was first used to map the PSF over 
the whole GMOS field of view, by determining the 
shape of unsaturated stars, which have been se- 
lected based on their FWHM. For the stars, we 
prevented im2shape to do any deconvolution by 
using a Dirac's delta function as the input PSF. 

^From the original sample of stars we kept only 
those that actually sample the local PSF. This has 
been done through a sigma-clipping process, where 
stars with too deviant ellipticities or position an- 
gles were removed. The remaining stars map a 
PSF which is nearly constant over the field, with 
average ellipticity of 0.056 (or 5.6%) and major 
axis oriented nearly East- West. 

The next step is to run IM2SHAPE for galaxies. 
As galaxy images have more complex shapes than 
stars, they were modeled by a sum of two two- 
dimensional elliptical Gaussians with the same 
center, ellipticity and position angle. For each 
galaxy, the input PSF was calculated using their 
five closest stars. 

3.2.2. Sample selection 

Once we have measured the galaxy shapes, we 
need to identify the background galaxies, which 
can be used as probes of the cluster shear field. 
These galaxies constitute what we call the weak- 
lensing sample. Since their redshift is unknown, 
they were selected by their magnitudes and colors 
(when available). We have included in the weak- 
lensing sample only galaxies fainter than r' = 23.5 
and with (r' - i') colors not closer than 0.2 mag 
to the cluster red sequence. These criteria are ad- 
equate because the number density radial profile 
of the weak-lensing sample is nearly flat and does 
not decrease with increasing radius, as one would 
expect if there was significant contamination by 
cluster members. 

In order to select a good quality weak-lensing 
sample, we kept only objects having errors on the 
tangentially projected ellipticities (with respect to 
the cluster center) smaller than 0.35. The re- 
sulting weak-lensing sample has 276 galaxies; the 
faintest has r' = 25.6 mag., and the average mag- 
nitude in this sample is r' = 24.8. 
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3.2.3. Mass distribution 

The information on the shape of these galaxies 
(and the corresponding errors) was used to feed 
the software LensEnt (Bridle et al. 1998; Mar- 
shall ct al. 2002) which, based on a maximum en- 
tropy method, reconstructed the projected mass 
density distribution of Abell 586. The resulting 
mass contours can be viewed in Figure 5. This 
map has been smoothed with a two-dimensional 
Gaussian with FWHM = 2.0', that maximizes the 
likelihood of the reconstructed mass density given 
the data. 

The mass distribution presented in this figure, 
showing a single mass clump associated with the 
BCG, is qualitatively very similar to the X-ray 
emission. This strongly suggests that the cen- 
tral region of Abell 586 is indeed in dynamical 
equilibrium. It is also worth noticing an exten- 
sion towards the SE direction of the cluster, which 
is present both in the X-ray and weak-lensing 
contours. Our map is also in qualitative agree- 
ment with the weak-lensing mass map produced 
by Dahle et al. (2002), despite of a small offset 
in the peak position which, however, is within the 
resolution of both maps. 




Right Ascension (J2000) 



Fig. 5. — GMOS image of the central region 
of Abell 586 superimposed with logarithmically 
spaced X-ray isophotes (continuous lines) and 
weak lensing reconstructed mass density (dashed 
lines). The X-ray point source near the southwest 
corner is the Seyfert 1 galaxy C171_3650. 



3.2-4- Radial shear profile 

Very often the shear profile is fitted by a singu- 
lar isothermal profile (SIS; e.g. Mellier 1999). This 
model is very convenient because it has a single 
parameter, the Einstein radius, 6e, and allows di- 
rect comparisons with results already obtained for 
this and other clusters. Besides, it offers a rough 
approximation when only the central part of the 
cluster can be accessed, as in the present case. In 
this model, 0e is directly related to the cluster 
velocity dispersion by the expression: 

e E = A7r^p, (i) 

c z 

where (3 is the ratio of Di s , the angular diameter 
distance between the lens and the source, over D s , 
the angular diameter distance to the source. 

We adopt a simple parametric model for the 
mass distribution, the singular isothermal profile 
(SIS). More complex models are not appropriate, 
because our data has a rather high noise level due 
to the small number of galaxies to probe the shear 
field. Additionally, this kind of profile is widely 
used in studies of this type, allowing a simple com- 
parison with other clusters. 

To estimate the average (3 of our sample we 
have used the catalog of photometric redshifts for 
the HDF of Fcrnandez-Soto, Lanzetta & Yahil 
(1999). From this catalog, we selected a sample 
with the same bright limit and average magni- 
tude as the present weak lensing sample. This 
catalog, however, does not have r' magnitudes; 
thus we have used I magnitudes, adopting the 
color (r'— I) = 1.55. This color is typical for a 
Sbc galaxy at z = 0.8 (Fukugita, Shimasaku & 
Ichikawa 1995). This procedure resulted in an av- 
erage value ((3) = 0.678, corresponding to an av- 
erage redshift of 0.78 ± 0.05. 

The best-fit shear profile (Figure 6) gives an 
Einstein radius of 30 ± 3", which translates into 
a = 1243 ± 58 km s -1 . The uncertainty of 58 
km s _1 is the statistical error of the fitting. It 
is important to mention that the determination of 
is the major source of systematic uncertainties 
on the derived velocity dispersion and thus on ab- 
solute mass determinations through weak-lensing. 
For instance, if we adopt the color of a typical 
Sab or Scd galaxy, instead of a Sbc, the resulting 
velocity dispersion changes about ±50 km s . 
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The value of the velocity dispersion we obtain 
here is significantly smaller than the one found by 
Dahle et al. (2002), 1680±ifg km s" 1 . ft is diffi- 
cult to figure out the reasons of this discrepancy. 
One possible reason is that Dahle et al. have used 
the approximation .9^7 instead of g = 7/(1 — k) 
used on the present work. Here, g is the distortion, 
which is directly related to background-galaxy el- 
lipticities, 7 is the gravitational shear, and k is the 
convergence or the projected mass density normal- 
ized with respect to a critical mass density: 



S 

si 



1 



4ttG Di (3 



(2) 



For strong lensing clusters like Abell 586, the 
assumption implicit in the approximation used by 
Dahle et al., i.e. n ~ 0, is too strong, mainly for 
the cluster central regions. However, adopting this 
approximation in our sample, we get a = 1353 
km s _1 . This valu is closer to the Dahle et al. 
result, but it is still significantly smaller. It is 
worth mentioning that the observations employed 
by them have similar depth to our observations. 
Dahle et al. observations were done with the 2.56- 
m NOT telescope with exposures of 5.4 ks in both 
V and / filters. After scaling and summing both 
exposures, the resulting imaging is about the same 
of the present 1.2 ks r' imaging with the 8.1-m 
Gcmini-N telescope. Both detectors have also a 
similar ficld-of-view. In terms of image quality, 
however, our data is better. Dahle et al. reported 
a seeing FWHM of 1.0 and 0.8 arcsec for their V 
and / images, respectively, whereas for our data 
this value is 0.7 arcsec. 

3.2.5. Total mass 

The mass profile can be directly measured from 
the shear data using the technique called aperture 
mass densitometry (AMD; Fahlman et al. 1994). 
This method is based on the single assumption 
that the surface mass density is circularly sym- 
metric. 

The AMD actually measures the difference be- 
tween the mass density inside a given radius and 
the ring between this radius and a maximum ra- 
dius of reference (r max ): 






] c 
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< 
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Fig. 6. — The shear profile for the cluster Abell 
586. The filled symbols correspond to average 
cllipticities of the faint background galaxies pro- 
jected tangcntially to the cluster center. Each 
point represents the average of at least 55 galax- 
ies of the weak-lensing sample. The open squares 
are the same, but for cllipticities projected in a 
direction 45° relative to the center. The solid line 
shows the best-fit SIS profile. The dashed line 
shows the Dahle et al. (2002) best SIS model (a 
= 1680 km s^ 1 ). 



where £ is defined by: 

C(ri,r max ) = k(r < r x ) - R{r x < r < r max ) 

= T^k^ ^r^r)d\nrM) 

where gt is the tangcntially projected distortion, 
r\ the radius in which the mass is measured, and 
R the mean convergence. 

We choose r max = 150 arcsec (437 kpc at z = 
0.17), which is the largest radius fully contained 
within the field of view. Therefore, the total mass 
inside 422 kpc (145") 4 estimated by AMD is (4.3± 
0.7) x 1O 14 M0. The mass profile can be seen, 
together with X-ray mass profiles, in Section 3.3.3. 

^From equation (3) it can be seen that the mass 
measured using the AMD method depends on S c 
and thus on (3. Therefore the same systematical 
uncertainties related to the poorly known rcdshift 
distribution of the background sources, as previ- 
ously discussed, applies here. 



Mamd(< r) 



2sr- 

7r r h r 



C(r, 



(3) 



4 largcst n so that there is at least 20 galaxies of the weak- 
lensing sample inside the annulus ri < r < r m ax' 
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3.2.6. Strong- Lensing Features 

The central region of Abell 586 (see Figure 7) 
shows several low surface brightness structures ori- 
ented tangentially to the cluster center, most of 
them with (g'-i') colors up to 0.5 mag bluer than 
cluster ellipticals with similar brightness. 

A prominent giant arc, already reported by 
Dahlc et al. (2002), can be seen in the Northwest 
direction, distant 22 arcsec of the BCG center. 
There are several high surface brightness galax- 
ies superimposed on this arc. At the opposite side 
of the central galaxy another giant arc, although 
fainter, can be appreciated in Figure 7 at 20 arcsec 
from the cluster center. Unfortunately, we could 
not determine the redshift of these arcs and, there- 
fore, could not confirm whether they are multiple 
images of the same source. 

We succeeded, however, to measure the red- 
shifts of two arclets in the vicinity of the BCG, at 
26.8 arcsec and 19.7 arcsec from the BCG center, 
whose orientations are nearly tangential to the di- 
rection of the cluster center. The spectra of these 
objects present emission and interstellar absorp- 
tion lines, both typical of late-type galaxies. Their 
measured rcdshifts are 1.43 and 0.61, respectively 
(cf. Figure 7). No other objects with colors simi- 
lar to those of these arclets were found, so no ad- 
ditional candidates for multiple images could be 
identified. 

In the absence of multiple gravitational images 
it is not possible to model the cluster potential 
using the position of the arcs (like in Kneib et al. 
1996). However, a rough estimate of the cluster 
mass inside the region enclosed by the arclets can 
be obtained assuming that they put a limit on 
the position of the Einstein radius. Under this 
assumption, we obtain, from the higher and lower 
redshift arclets, a equal to 1056 km s _1 and 998 
km s^ 1 , respectively. 

This cluster shows another strong-lensing-like 
feature that deserves be mentioned. In Figure 8 
we present a close-up of the spiral galaxy and clus- 
ter member C281_3813, that is at 1.7' from the 
cluster center. What is remarkable in this figure 
is the presence of an arc-like object (indicated by 
an arrow). If this object has indeed a gravita- 
tional origin, it would be an uncommon example 
of strong-lensing by a late-type galaxy. 
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Fig. 7. — Optical image of central region of Abell 
586. Objects with redshifts measured by us inside 
this field are marked. The dashed lines surround 
the gravitational arcs. This image shows a square 
region of 1.2' on a side. 
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Fig. 8. — Image of the galaxy C281_3813, showing 
an arc-like object. 



3.3. X-ray 

3.3.1. X-ray brightness profile 

The X-ray brightness profile of Abell 586 was 
obtained with the STSDAS/IRAF task ellipse. 
We have used the exposure map corrected image 
in the [0.3-8.0 keV] band, binned by a factor 4 (one 
X-ray image pixel has 2 arcsec). We have masked 
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the CCDs gaps and source points. The brightness 
profile, shown in Figure 9, could be measured up to 
500 arcsec (1.46 h? Mpc) from the cluster center. 
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Fig. 9. — X-ray brightness profile. The full line 
shows the best fit /3-model with = 0.518 ±0.006 
and R c = 23.1 ± 0.6 arcsec (67 ± 2/i7 1 kpc); the 
dotted line shows the /3-model plus a constant 
background brightness (which was actually fitted) . 
The vertical dotted-dashed line indicates the virial 
radius (cf. Section 3.3.3). 

Using the /3-model (Cavaliere & Fusco-Femiano 
1976) to describe the surface brightness radial pro- 
file, 



S x (i?) = S n 



-3/3/2 



(5) 



a least-squares fitting gives (3 = 0.518 ± 0.006 and 
R c = 23.1 ± 0.6 arcsec (67 ± 2/if 1 kpc). If we 
assume that the gas is approximately isothermal 
and distributed with spherical symmetry, there is 
a simple relation between the brightness profile 
and the gas number density, n(r), i.e., 



n(r) = n 



1+1- 



21 -3/3+1/2 



(6) 



where R c — r c (capital indicates projected 2D co- 
ordinates, lower case indicates 3D coordinates). 

In order to estimate the central electronic den- 
sity, no, we integrate the bremsstrahlung emis- 
sivity along the line-of-sight, in the central re- 
gion, and compare with the flux obtained by spec- 
tral fitting of the same region. We obtain thus 



3.3.2. Temperature profile 

We have computed the radial temperature pro- 
file using concentric circular annuli. For each an- 
nulus, defined by approximately the same number 
of counts (2000 counts, background corrected), a 
spectrum was extracted and fitted following the 
method described above, except that the hydro- 
gen column density and metallicity were kept fixed 
at the mean best- fit value found inside 1.4 arcmin 
(i.e., N H = 9.lti j x 10 20 cm~ 2 and the metallicity 
Z = 0.5 Zq). Figure 10 shows the temperature 
profile. 

R [ft^kpc] 
400 600 800 
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ACIS-I 




400 



Fig. 10. — Temperature profile. The error bars are 
la confidence level. The horizontal dotted-dashed 
line is the mean temperature inside 1.4 arcmin, 
kT = 7.8lJ'g keV, the horizontal dotted lines cor- 
respond to 90% confidence level of the mean tem- 
perature. The full line shows the best-fit of the 
poly tropic model (see text). 



(18.4 ± 1.5) x 10" 3 cm- 



Since the temperature profile shows a clear gra- 
dient, we have tried a simple parametric form of 
the ICM temperature using a polytropic equa- 
tion of state. It is not clear whether a polytropic 
model describes well the ICM. Indeed, Irwin et 
al. (1999) and Dc Grandi & Molendi (2002) ar- 
gue that this is not a good description of the gas 
temperature profile in clusters. However, hydro- 
dynamic simulations (Suto et al. 1998), theoret- 
ical models (e.g., Komatsu & Seljak 2001; Dos 
Santos & Dore 2002) and some observations (e.g., 
Markevitch et al. 1999; Finoguenov et al. 2001a) 
suggest that the gas may be described by a poly- 
tropic model, with polytropic index 7 roughly in 
the range 1.1 < 7 < 1.2. 
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Therefore, we have fitted a polytropic temper- 
ature profile: 



TO) = To 



-, -3/3( 7 -l)/2 



(7) 



where r c and (3 are the values obtained with the 
/3-model fitting of the brightness surface profile, 
and To is the central temperature. 

A standard least-square fit of Eq. (7) with only 
two free parameters, To and 7, yields a rather good 
fit: T = 8.99 ± 0.34 keV and 7 = 1.10 ± 0.03, cf. 
Figure 10. This value agrees with those found by 
Finoguenov ct al. (2001b) and, being well below 
5/3 (the ideal gas value), suggests that the gas 
may be in adiabatic equilibrium (see, eg., Sarazin 
1988, Sect. 5.2). 

We note that this cluster does not present any 
sign of cooling in the very central part, at R w 
70 /itq 1 kpc, the smallest radius where we can ex- 
tract a meaningful spectrum and measure the tem- 
perature. Either we lack the resolution to detect 
an eventual drop in temperature or the intraclus- 
ter gas is not cooling. Since the central cooling- 
time is roughly 



tcooi ~ 5.8 x 10 



T i/2 

9fkoV 

n 3 



10 years , 



(8) 



then, if the gas is indeed not cooling, something 
must be heating it (as it was already realized for 
some clusters, e.g., Peterson et al. 2001). Heating 
by cluster merging seems improbable, given the 
apparent spherical symmetry of the X-ray emis- 
sion. Other possibilities, like heating by AGN, 
thermal conduction, etc., may be playing a role in 
the energy budget of this cluster (e.g., Markevitch 
et al. 2003). 

However, we may be simply not detecting an 
eventual drop in temperature because we lack the 
resolution. Using a sample of 20 clusters, Kaas- 
tra ct al. (2004) show that the radius (r c in their 
paper) where the temperature drops in cooling- 
flow clusters is, with 2 exceptions, smaller than 
70/ifo 1 kpc. 

3.3.3. Gas and total masses 

We estimate the gas mass simply by integrat- 
ing the density given by Eq. (6), assuming spher- 
ical symmetry which, in this case, seems a good 



approximation (cf. Figure 5). The integrated gas 
mass is shown in the top panel of Figure 1 1 . 
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Fig. 11. — Top: X-ray dynamical mass profile 
assuming isothermal profile (full line), polytropic 
profile (dashed), weak-lensing mass profile using 
the AMD method (small dashed) and strong- 
lensing mass estimates (triangles). The shaded ar- 
eas represents l-a uncertainties. The arrow on the 
top of the panel marks the radius that contains 
all the galaxies used for the dynamical analysis. 
Bottom: Baryon fraction for isothermal (full line) 
and polytropic (dashed) temperature profiles. The 
vertical line indicates the virial radius. 

Having the temperature profile, we also com- 
pute the "X-ray" dynamical mass (i.e., the dynam- 
ical mass determined from an X-ray observation, 
not to be confused with the mass of the X-ray 
emitting gas). Supposing hydrostatic equilibrium 
and polytropic temperature profile, the dynamical 
mass is given by: 



M(r) 



3kT f3jr c ( r 



-l-§(7-l)0 



(9) 



If 7 = 1, we have the usual isothermal mass pro- 
file. Here, fj,m p is the mean mass per particle 
(fx = 0.6 for a fully ionized gas with primordial 
He abundance) and m p is the proton mass. The 
top panel of Figure 1 1 shows the dynamical mass 
profile estimated both with an isothermal (7 = 1) 
and polytropic temperature profile (7 = 1.1). 
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The gas mass at R = 400" (1.16 /i^Mpc) is 
0.43 x 10 14 M©. Depending on the assumed tem- 
perature profile, A/dyn = 5.3xl0 14 M Q (isothermal 
case) and Md yn = 4.25 x 10 14 M Q (polytropic case) 
within the same radius. 

Virial radius and baryon fraction 

We can now compute the radius corresponding to 
the density ratio 6 = p(rs)/p c (z). For 5 = 200, 
we have the usual r2oo; for the SCDM and ACMD 
cosmological models, 5 = 180 and 340 correspond, 
respectively, to the virial radius (cf. Lacey & Cole 
1993; Bryan & Norman 1998). 

Taking into account the polytropic temperature 
profile, we use the formula in the appendix of Lima 
Neto et al. (2003) to compute the various r$. Ta- 
ble 4 shows the results, where each column corre- 
sponds to a density ratio and the lines with "poly" 
and "isoth" refer to the polytropic and isothermal 
temperature profiles, respectively. Then, assum- 
ing a polytropic profile and the ACDM scenario, 
the virial radius is r v j r = 1.2&hjQ Mpc, which cor- 
responds to R = 430". 

Using the dynamical and gas mass profiles, we 
can derive the baryon fraction profile, f(r) = 
Afbaryon/Aftotai, assuming that the bulk of the 
baryons are in the ICM (i.e., Mbaryon ~ M gas ). 
The galactic contribution to the baryon mass is 
taken into account following White et al. (1993) 
and Allen et al. (2002): 

M gas (l + 0.16/4 5 ), (10) 

where M gas is the intracluster gas mass. 

At R = 400 arcsec, the baryon fraction is still 
rising. Its value at this point depends on the as- 
sumed temperature profile: for the isothermal and 
polytropic cases we have / = 0.10 and 0.12, re- 
spectively. The baryon fraction determined here 
agrees with the usual values found in rich clusters 
(e.g. White & Fabian 1995; Allen et al. 2002). 

At the virial radius, defined above, we have 
M gas = 0.48 x 10 14 M Q , M dyn = 4.53 x 10 14 M Q , 
and / = 0.12 (assuming a polytropic temperature 
profile). 

4. Discussion 

The mass determinations resulting from the ap- 
plication of four distinct techniques, based on dif- 



ferent types of data, to the cluster Abell 586 are 
compared in this section discussed within the con- 
text of eventual deviations from a relaxed state. 

In Table 5 are listed the velocity dispersions ei- 
ther measured or deduced in this paper through 
strong-lensing, X-ray, redshift survey, and weak- 
lensing methods, together with results of other au- 
thors. 

Except for the velocity dispersion determined 
by Bottini (2001) and Dahlc et al. (2002), all es- 
timations agree well within at least a 68% confi- 
dence level. This result is strengthened when we 
compare the mass profiles provided by each tech- 
nique (see Figure 11), although systematic differ- 
ences within 2<T between the (model dependent) 
measurements are noticed. 

Recently, Cypriano et al. (2004) compared mass 
estimates obtained through different techniques 
for a sample of 24 X-ray luminous clusters with 
z < 0.3 and with homogeneous weak-lensing deter- 
minations. Adopting the criteria that the agree- 
ment, or disagreement, between dynamical (ve- 
locity dispersion and X-ray) and non-dynamical 
(weak-lensing) mass estimators in a particular 
cluster is indicative that the cluster is relaxed, or 
not, we have found that clusters with ICM tem- 
peratures above ~8.0 keV show strong evidences 
of dynamical activity, while cooler clusters tend to 
be nearly relaxed. 

Our study of this cluster suggests that Abell 
586 is a well relaxed object that has not experi- 
enced a major merger in the last few Gyr; notice 
that its ICM temperature is just below the upper 
~8.0 keV limit found for quasi-relaxed systems, 
so that this indicator may have a valid predictive 
character. 

Gravitational arcs are often found in clusters 
that are not relaxed. This behavior can be un- 
derstood not only because the more massive clus- 
ters are in general young structures, but also be- 
cause the presence of substructures and other fea- 
tures associated with dynamical activity, enhance 
the strong-lensing cross-section, as shown, for ex- 
ample, in ray-tracing simulations by Bartelmann, 
Steinmetz & Weiss (1995). 

Indeed, a recent study by Smith et al. (2004) 
of a sample of 10 strong-lensing clusters, selected 
by X-ray luminosities, using HST and high quality 
X-ray data, concluded that only 30% of them can 
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Table 4 

Radius at different density contrast. 



8: 180 


200 


340 500 2500 


Model 


r s : 1.69 


1.61 


1.26 1.05 0.49 


poly 


arcsec 1.97 


1.83 


1.40 1.16 0.51 


isoth 


Obs.: 


r s is 


in units of hj Mpc. 





Table 5 

Cluster velocity dispersion derived with different methods. 



(1) 


(2) 


(3) 


a [km s : ] 


Method 


Notes 


998-1056 


strong-lensing 


9 E 19.7"-26.8", sources z = 0.61 and 1.43 


1015 ±500 


X-ray luminosity a 


L x = 18 x 10 44 erg s" 1 


1174 ± 130 


X-ray temperature 3. 


kT = 7.8 keV 


1050 ± 350 


X-ray 11 


Allen (2000) 


1161 ± 196 


velocity distribution 


31 galaxies 


313 ± 70 


velocity distribution 


7 galaxies Bottini ( 2001 ) 


1243 ± 58 


weak lensing 




1680 ± 170 


weak lensing 


Dahle et al. (2002) 



a For X-ray data we have used the empirical relations Lx~o~ and Tx~o~ from Xue 
& Wu (2000). 



Note. — Italic entries are quoted results from the literature. 



be classified as truly relaxed clusters. Actually, 
the Cypriano et al. criteria successfully predicts 
the dynamical state of 80% of this sample. 

Given that cluster inner mass profiles are often 
determined by strong- lcnsing analysis (e.g., Tyson, 
Kochanski & dell' Antonio 1998; Smith et al. 2001; 
Sand, Treu & Ellis 2002; Kneib et al. 2003; Sand 
et al. 2004), and possibly the majority of these 
systems are non-relaxed, then, the disagreement 
between these observed profiles and the theoreti- 
cal dark matter profiles, derived from relaxed ha- 
los found in numerical simulations (e.g. Navarro, 
Frenk & White 1997; Ghigna et al. 1998; El-Zant 
et al. 2004), might be due to the fact that the 
physical state of the observed and modelled sys- 
tems are not consistent. 

The identification and detailed mass recon- 
struction of a representative sample of clusters 
like Abell 586, using several techniques like grav- 
itational lcnsing, X-ray emission and galaxy ve- 
locities, is a promising way towards a better un- 
derstanding of the behavior of baryonic and dark 
matter components in the center of galaxy clus- 
ters. 

5. Conclusion 

Using optical data taken with the 8-m Gemini- 
North telescope and available Chandra X-ray 
data for the A586 galaxy cluster, we have de- 
rived its mass distribution and content using a 
multi-technique analysis. Our main results can be 
summarized as follows: 

1. Radial velocity measurements for 31 cluster 
galaxies resulted on a systemic redshift of 
(z) = 0.1708 ± 0.0001 and a velocity disper- 
sion a = 1161 ± 196 km s _1 . 

2. We detected weak gravitational shear, whose 
best fit through an isothermal mass profile 
gives a = 1243 ± 58 km s _1 . 

3. We identified a system of gravitational arcs 
and determined the redshifts for two arclets 
(z = 0.61 and 1.43) belonging to this system. 

4. We determined the mass distribution in the 
central region of the cluster through two 
techniques: wcak-lensing and X-ray emis- 
sion; they are found to be very similar and 
having almost circular geometry 



5. The ICM gas is distributed very smoothly; it 
has a mean temperature of 7.8^' g keV and 
a mean metallicity of O^I^q'^Zq, both val- 
ues being slightly higher than the averages 
reported in the literature for rich clusters. 

6. The gas temperature profile is well described 
by a poly tropic model with 7 = 1.1. 

7. The cluster virial radius is approximately 
l.ShjQ Mpc, and the gas and dynami- 
cal mass within this radius are M gas = 
0.48 x 10 14 M© and M dyn = 4.53 x 1O 14 M ; 
the baryon fraction at the same radius is 
/ = 0.12, assuming a polytropic tempera- 
ture profile. 

The ensemble of our observational results, de- 
rived with different techniques and wavelcnght 
ranges, indicate consistently that A586 is a mas- 
sive cluster, characterized by a velocity dispersion 
that is in the range 1000-1250 km s _1 . Several 
pieces of evidence suggest that this cluster is dy- 
namically well relaxed, namely: i) The near circu- 
lar mass and X-ray luminosity distribution, both 
concentric with the BCG; ii) The agreement, un- 
certainties taken into account, between dynami- 
cal (X-ray, galaxies velocity dispersion) and non- 
dynamical (gravitational lcnsing) mass estimators; 
in) An ICM temperature profile well described by 
a polytrope with index 7 = 1.1. 

As a final remark, it is interesting to note that 
Abell 586 is found to follow the Cypriano et al. 
(2004) criteria to diagnose the dynamical state of 
luminous X-ray clusters. Its ICM temperature is 
just below the upper ~8.0 keV limit claimed by 
Cypriano et al. for quasi-relaxed systems. 
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